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We present a first-principles derivation of the Markovian semi-group master equation without invoking the 
rotating wave approximation (RWA). Instead we use a time coarse-graining approach which leaves us with a free 
timescale parameter, which we can optimize. Comparing this approach to the standard RWA-based Markovian 
master equation, we find that significantly better agreement is possible using the coarse-graining approach, for 
a three-level model coupled to a bath of oscillators, whose exact dynamics we can solve for at zero temperature. 
The model has the important feature that the RWA has a non-trivial effect on the dynamics of the populations. 
We show that the two different master equations can exhibit strong qualitative differences for the population 
of the energy eigenstates even for such a simple model. The RWA-based master equation misses an important 
feature which the coarse-graining based scheme does not. By optimizing the coarse-graining timescale the latter 
scheme can be made to approach the exact solution much more closely than the RWA-based master equation. 



I. INTRODUCTION 



The reduced dynamics of a quantum system coupled to another quantum system (the environment, or bath) can always be 
described by a completely positive trace preserving (CPTP) map provided the initial state of the system and the bath is factorized 
[1 1, or classically correlated [2 3 |. The map can always be formally written in the Rraus operator sum representation (OSR): 

p(t) = Y / Mt) P (0)A\(t), (1) 

i 

but calculating the Kraus operators Ai explicitly is usually too difficult, making this form of the map unusable for most practical 
purposes. It is often more convenient to work with a differential equation that describes the evolution. Many such "master 
equations" have been derived, either from first principles or phenomenologically, and have been hugely successful in describing 
a wide range of phenomena fl]|4r|7]]. However, master equations in general do not necessarily generate CPTP dynamics, and can 
even violate the positivity of the density matrix for certain initial states. Lindblad showed that a Markovian semi-group master 
equation (SME) that is guaranteed to generate a CPTP map must be of the form: 



ip(t) = -T[H, P (t)] + |^ 7a/3 ([i Q ,p(i)4] + \L a p(t),Ll 



dt ry ' h 1 ' rwj 2 



(2) 



where H is Hermitian and the coefficients jap form a positive matrix. All first-principles derivations of the master equation that 
attempt to end in this form require what is often called the "rotating wave approximation" (RWA), also referred to as the "secular 
approximation". This approximation requires dropping terms which occur on timescales of order l/(ui — u/) <C tr, where uj 
and to' are system energy eigenstate differences and tji ~ 1 /g is the typical relaxation time in the interaction picture, where g 
is the system-bath interaction strength. The RWA is notoriously hard to justify rigorously and is well known to sometimes give 
wrong results, e.g., Refs. 18 UT21 . Nevertheless, it has been used extensively in the literature, and so it is tempting to formulate a 
general alternative approach to the first principles derivation of Markovian master equation. 

In Refs. lfl3l[T4l . a (phenomenological) method was proposed to derive a Markovian SME from any CP map without requiring 
the RWA. The procedure involves averaging the dynamics over a coarse-graining timescale. Therefore, we can derive two 
SMEs, one using the first-principles derivation with the RWA (which we refer to as the RWA-SME) and the second using the 
phenomenological coarse-graining idea of Refs. lfT3l[T4ll (which we refer to as the CG-SME). In order to properly compare the 
effect of the RWA, we require a model to which an exact solution can be found, and in which the RWA acts non-trivially. We 
construct a three-level system coupled to an infinite harmonic oscillator bath, the dynamics of which we can solve exactly at zero 
temperature. Using this model, we compare the results of the two SMEs to the exact solution. We find that even in this simple 
example, the dynamics of the RWA-SME misses an important qualitative feature of the populations of the exact solution, which 
the CG-SME does not. Moreover, no attempt has been made so far to optimize the coarse-graining timescale, and we show that 
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the CG-SME can be fine-tuned to better match the exact solution via such an optimization. 

The structure of this paper is as follows. In Section [II] we first review the coarse-graining approach of Refs. Ifl3l [141. then 
derive a coarse-grained SME from first principles starting from the cumulant expansion approach of Alicki |[T5l . In Section[HI]we 
describe an exactly solvable 2-level model where the RWA-SME and CG-SME do not make substantially different predictions. 
Motivated by this observation, we construct and solve a 3-level model in Section IV where the two master equations do make 
substantially different predictions, and derive the parameters of the two master equations in the context of the model in Section|V] 
We then compare the predictions to the exact solution in Section |VI[ where we show that by optimizing the coarse-graining 
timescale we can get better agreement with the exact solution by using the CG-SME. In Section VII we briefly discuss the 



question of how to optimally choose the coarse-graining timescale in a more general setting, and conclude in Section VIII 
Various technical details are presented in the Appendix. 



II. SEMIGROUP MASTER EQUATIONS VIA COARSE GRAINING 
A. From CP maps to the SME 

In this subsection we briefly review the derivation of the SME |[T3l [141 . to motivate the idea of time coarse-graining. We 
consider the dynamics of a system S coupled to a bath B, together forming a closed system. The closed system evolves unitarily 
under the Hamiltonian 

H = H S + H B +H I , (3) 

where Hs is the system Hamiltonain, Hb is the bath Hamiltonian, and Hi is the interaction Hamiltonian. Let us work in the 
interaction picture (we will denote interaction picture operators with a tilde (" " "), such that the total evolution operator is given 
by: 



U(t) = T+ exp 



H I {s)ds 



o 



where 7+ denotes the time-ordering operator. The system evolution is given by: 

ps(t) = Tv B \u(t) Ps (0) ® p B (0)U\t) 



(4) 



(5) 



where we have assumed that at time t = the closed system density matrix is factorized. By writing the spectral decomposition 
in the basis of the bath density matrix eigenstates, 

PB (0) = ^A»( M |, A m >0, 5>„ = 1, (6) 
n n 

we can write the evolution equation in the operator-sum representation (OSR): 



PS 



(t)=J2Mt)ps(0)Al(t), (7) 



where the Kraus operators {A{\ are given by 

A i (t) = yfc{lAU(t)\v), i=(n,v). (8) 

The Kraus operators have the property: 

= ^2x v {u\W\^\U\y) = £A„ = 1 . (9) 

It is useful to expand the Kraus operators in terms of a fixed (time independent) operator basis {K a }^ =0 with Kq = 1 and 
Nb the dimension of the bath Hilbert space, where M = N B — 1, such that: 

M 

A i (t) = ^2b ia {t)K a , (10) 

a=0 
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to give: 



M 



Ps(t) = E X*p(t)K aPs (0)Kl 

a,/3=0 



(ID 



with Xap(t) = J2i=i b i*(t)b* if3 (t). In particular, X*p a (t) = Xap(t) and \ aa = Ei \ka\ > so X is Hermitian with positive 
diagonal elements. The property in Eq. |9]) now gives 



M 



i a, f3=0 



(12) 



By multiplying this expression from the left and right by p(0)/2 and adding the two expressions, we get the following relation: 



PS 



M M M 

(O)=Xoo(t)ps(O) + -J2x a o(t){K a ,p s (O)}+-J2xo0(t){ps(O),Kl} + - ]T X *p(t) {k\k*, psQ)} . (13) 



I3=i 



a,/3=l 



Taking the time derivative of Eq. ( fTT) , we are able to write an evolution equation: 



E WW K a ,p s (0)Kl + K aPs (0),K i 



a, 13=1 



where we have used the time derivative of Eq. ( fT3| l to replace terms with Xoo(0 an d defined 



S(*) = ?E(Xa0(*)#a-X0a(i)#i) 



(14) 



(15) 



Although this form of the evolution equation has a striking resemblance to the SME in Eq. pi, it determines ps{t) in terms 
of ps(0), i.e. at a special time t = where the state is factorized into a product state. To write an evolution equation such 
that the equation depends on ps(t), we must make some approximations. Let us define three different time scales: (i) a time 
r c corresponding to the inverse of the high-frequency cut-off ui c in the bath spectral density, (ii) a system time scale t$ cor- 
responding to the characteristic time for significant changes in the system density matrix in the interaction picture, (iii) and a 
coarse-graining time scale At corresponding to the time-scale after which the bath has effectively "reset." Let us assume the 
following relationship holds: 



T c < At < T S . 

Let us define p n — p(nAt), such that discretizing Eq. ([14} in terms of At gives: 



Pi - Pa 
At 



(S),Po 



aP) 



K a ,p K 



K a p ,K 



where 



(X) 



At 



At 



X(s)ds 



(16) 



(17) 



(18) 



In principle, Eq. ( p~4] > cannot be used to determine the next discretization step relating p2 in terms of p\ . However, if we assume 
that the bath resets in the time A< such that the bath interacts with the system in exactly the same manner at each time step At, 
then we can simply proceed with our discretization. This amounts to performing a Markovian approximation, and we are led to 
the form of the SME: 



d_ 
dt' 



M 

Ps(t) - -\ [(S),p s (t)] + 2 E <*«/>) + [K*Ps(t),K P ]) . 



(19) 



a,P=l 
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It was shown in Ref. [14] that the coefficients (Xap) form a positive matrix, which means that Eq. ( fT9| ) is in Lindblad form 
and hence generates completely positive dynamics. At no point was the RWA made in this derivation. 

Note further that the coarse-graining timescale At enters Eq. ( fT9| i in a fundamental manner via the coarse-grained averages 
(S) and (Xa/3) ■ Clearly, Eq. ( [To} leaves room for optimization of At, a point we shall return to in our discussion below. 



B. From the cumulant expansion to the SME 

We now present a semigroup master equation that can be derived from a first-principles theory, again without invoking the 
RWA, using coarse-graining instead. The derivation is somewhat lengthy, so we present it in four parts. We follow Alicki's work 
lfT51 for the first two parts. 



1. Cumulant expansion 

We first review the cumulant expansion approach of Ref. ITSl . Let us consider a Hamiltonian: 

H = H S + H B + XHi 

with 

Hi = A ® B 



(20) 



(21) 



where A is a Hermitian system operator and B is a Hermitian bath operator. We have restricted ourself to a single term to 
simplify the notation, but the more general case with multiple terms follows in an analogous fashion. Let us assume that: 



where 



(B) B =0, p(0) = ps {0)®p B 
(X) b =Tt[ Pb X], 



and where ps is the Gibbs state for the bath: 



-0H E 



Pb 



Z B ' 

with Zb = Tr[exp(— /3Hb)] the partition function. Define: 

H Q (t) = H s {t) ®I B +I S ®H B 
U (t,0) = Texp (-i H (s)ds 
p(t) = ul(t,0)pU (t,0), 
where p(t) is the state in the Hq -interaction picture. In the interaction picture we have: 

d 



dt' 



P {t) = -i Ul\HiU ,p 



XHi(t),p 



with: 



Hi(t) = ul(t)AU s (t) ® U B (t)BU B (t) , 



where Us(t) is the unitary generated by H$- The formal solution is given by: 



p(t)=p(0)-i [ ds\xHi(s),p(s)] = p(0)-iX [ ds ifj (*)./)((>) 



(22) 
(23) 

(24) 

(25a) 
(25b) 
(25c) 

(26) 
(27) 



A 2 / ds ds' 
Jo Jo 



Hi(s), Hi(s'),p(0) 



(28) 
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We are interested in the reduced density matrix: 



ps(t) =Tr B (p(t)) = W t ps(P) 



(29) 



The cumulant expansion is given by 



W t = exp (^>2 A "^ (n) = 1 + XK(1) (*) + A " (^ (2) (*) + \ ( K(1) (*)) 



and we solve for by matching powers of A with Eq. (28 I. We get: 



0(A 3 



(30) 



KW(t)pg(Q) = -i f dsTr B ([HP(s),p(0) 



0, 



(31) 



where we have used (B) B = (which can always be done by shifting the operator B without loss of generality), and have used 
the shorthand: 



A(t) = U\AU S , B(t) = U B BU B , H { j L> (t) = A(t)® B(t) . 



(32) 



The next order in A gives: 



K^(t)p s (0) = -J o ds J ds' Tr B [[H^(s), [H^(s'),p(0) 
(Note that our expression differs by a factor of 1/2 in the second term relative to Eq. (4.10) in Ref. [15 1.) We write: 



frW, 



Tr, 



&P(s), H^(s'),p(0) ) = (A(s)A(s')p s (0) - A(s')p s (0)A(s))Tr B (B(s)B(s')p B ) + h.c. 



(33) 



(34a) 



(A^s)A(s')p s (0) - A(s')p s (0)A^s)) Tr B (B^s)B(s')p B ) + h.c. (34b) 



A simple norm estimate shows that ||i<"( n ) | = 0[(gt) n ], where g — \\Hi\\. Therefore the terms of third order and above can be 
neglected provided Xgt -C 1. This amounts to making the Born approximation. 



2. The second order cumulant 



We restrict ourselves to the case of time-independent Hg, which was also considered in Ref. |[T5"1 . Here we give a much more 
detailed derivation. Let 



A(t) = ul(t)AU s (t) =Y,e- l{£b - £a)t \e a )(e a \A\e b )(e b \ =J2^e- 



a . b 



such that 



A(s') = ^2A u e 

UJ 

This then gives the following map from time to t: 



K^(t)ps(0) 



]T (AlA u p s (0) - A u p s (0)Al) B uu ,(t) + h. 



where 



B uu ,{t)= / ds d s 'e l ^' s -^B( S , S ') , 
Jq Jo 

B(s,s') = (B\s)B(s>)) = (B(s)B(s , )}=B(s , ,sy 



(35) 



(36) 



(37) 

(38a) 
(38b) 
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We show in Appendix [A] how Eq. ( |37| i can be rewritten as 

K^(t)p s (0) = -i [S(t),p s (0)} + ]T A w p s (0)Al - - {aI,A u ,ps(0)} 

UJ,Lk)' 

where 

b uu ,(t) = f ds [ ds'e i ^' s -^B(s,s') = b* u , bJ (t). 
Jo Jo 

The "Lamb shift" is 

S{t) = Y J S^{t)Al,A w , 

where 

S^{t) = - l -J^ ds J ds' (e^' s -^B{ SlS ')-e- l ^ s -"'^B{s',s)) . 
Note that (S uu ')* = S^u, so that = S, as required for the interpretation of S as a Hamiltonian. 



(39) 



(40) 



(41) 



(42) 



3. Complete positivity 



Clearly, the dissipative (second) term on the RHS of Eq. ( |39| l appears to be in Lindblad form, but we must still prove the 
positivity of the matrix b(t) (this was not done in Ref. |15|). To this end we again expand the bath density matrix in its 
eigenbasis [Eq. |6|], and use this to write the correlation function B(s, s') = (£^(s)i3(s')) b explicitly. Let v be some arbitrary 
vector; then positivity amounts to showing that w6(<)t^ > for all v. Indeed: 

ifc(t)^ = £i;Jw(*K,= I ds [ ds' ^( Ua) e- W )^K/e-^'TTr[X]A M |M}<M|St (s)B(s')] (43a) 
, „ ,' Jo Jo , , , ,, ., 

= J2^ OA F\t)F(t) = ]T \p) || 2 > 0, (43b) 



where F(t) = g dsB(s) ^ v^e^ . 

Note that we have arrived at Eq. (|39]l without making the Markov approximation. Therefore, our (truncated) non-Markovian 
CP map is given by: 

ps(t) = e x2K(2) ^ Ps (0) (44) 
The complete positivity of this expression follows from the Lindblad theorem, as K^ 2 '(t) is in Lindblad form for all t > 0. 



4. SMEfrom the cumulant expansion and coarse-graining 

Finally, we show how to obtain the SME from the results above (this was not done in Ref. [ 15 1 either). Consider the CP map 
Eq. |44} and write it in the form of an evolution equation. Expanding the exponential to second order in A and taking a time 
derivative, we have: 



j t p s (t) = [x 2 s,p s (o)] + ]T \ 2 b^, k,Ps(o)4< - \ {aIa u ,p s (o)} 



Written in this form, the equation looks identical in form to Eq. ([14}, so we can make the identification: 



(45) 



X'b„ 



K a ^A u , K\ -> Al 



(46) 
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As discussed in section 
In particular, we have: 



II A the path to a SME is simply to take the time average of S and \ over the coarse-graining time At. 

\ 2 {S) = ~S(At)=Hl s , (47a) 
A 2 

X 2 (L,w) = & w , w >(Ai) = . (47b) 

The generator has the free parameter At that has to be fixed by some means. 

It is useful at this point to compare the SME derived above with that typically derived using the RWA in the time-independent 
case. For the latter case, which we refer to as the "RWA-SME", the SME takes the form (TJ: 



dt 
with 



d p s (t) = -t[H LS ,p s (t)} +$>M [AupsWAl - \ {AlA U) p s (t)}) (48, 



/oo 
dse ws B(s, 0) , H LS = S(oj)AIA u , 
- OO 



(49a) 



S(u) = / du/'y(u/)V 7 • (49b) 



(50) 



In the former case, which we refer to as the "CG-SME", the system operators A u are the same 

j t ~ Ps (t) = -i[Hl s ,p s (t)] + ^ 7ww /(At) (A u p s (t)Al - X - [AlA^psit)} 

but the rates 7 keeps a dependence on two different eigenfrequencies lj and uj f : 

7^'(At) = ~bu,u,'(At) , Hl s (At) - J2 S^.{At)Al,A u , (51a) 

UJ.LJ 1 

S ^'( At ) = -^ t J Q ds j o ds ' (e^' s ^B{ SlS ')- e-i^'-^BW ; S )) . (51b) 

We show in Appendix [B] that the RWA result is recovered by taking the At — > oo limit of the above result. This is not 
surprising, since the RWA (as it is performed in the standard first-principles derivation of the Lindblad SME) is equivalent to 
averaging the time-dependent ME operator over several periods of the oscillating terms. In Ref. [ 14| it was argued via a spin- 
boson example that ultimately the CG-SME generator should not depend on the averaging time At, and this is true if At is well 
beyond the bath correlation time tq. This corresponds to the RWA limit as mentioned above (i.e., At — > oo), at least for the 
SME generator that results from the CP-map derived from the cumulant expansion. As we desire to derive a master equation 
that accounts for dynamics slightly beyond the RWA regime, we do not take the RWA limit and keep the dependence on the 
coarse-graining time scale. 

Our main focus will be a comparison of the predictions of Eqs. ( |48| ) (RWA-SME) and ( |50| > (CG-SME), both to each other, and 
to the exact solution of a model we describe next. 



III. AN EXACTLY SOLVABLE MODEL WITH TRIVIAL ROTATING WAVE APPROXIMATION 



Let us consider the problem of a single 2-level system coupled to a bath of harmonic oscillators: 

H S = ■■ R B = J2 Wfe5 I &fe ' XHSB = X az ( b k + h ) ' (52) 

k k 

whose exact solution is known [ 1 ]. We wish to compare the exact solution with the solution found using the RWA-SME and the 
CG-SME. In all three cases, the diagonal elements of the density matrix are independent of time, and the off -diagonal component 
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takes the form: 

Pi 2 W = e^ ot e r( *Vi 2 (0). (53) 
For the three cases, the dephasing rate is given by: 

IW*) = -4A 2 f ^JM coth (^) (^^) (54) 
2 t f 00 , T . . , //3uA /l-cos(wAi) 



r CG (t, At) = -4A^— / dw J '(w) coth ^— £ i (55) 



r RWA (t) = -47rA 2 i lim ( - J(t j ) (56) 



where r exact is the exact result, Tcg is the CG-SME result, and Trwa is the RWA-SME result (details of the derivation are given 
in Appendix |C|). Formally, Tcc{t, t) — r exact (t). However, recall that At is supposed to be a fixed timescale, so setting At = t 
isn't in the spirit of the CG approximation. Both the CG- and RWA-SME yield a linear time-dependence for the dephasing, 
which is to be expected because of the Markovian approximation. Furthermore, the CG-SME reduces to the RWA-SME in the 
limit of At —> oo (also shown in Appendix O. Since the SMEs have the exact same behavior for the populations of the two 
energy eigenstates as the exact solution, we refer to this as a trivial case of the RWA. 



IV. EXACTLY SOLVABLE MODEL WITH NONTRIVIAL ROTATING WAVE APPROXIMATION 

We have so far presented a derivation of a SME without the RWA. In order to investigate the effect of taking/not-taking the 
RWA, we require an exactly solvable model that has more than one bath-induced transition. To this end, we consider a three 
level system (with energy eigenstates denoted by |0), |1), and |2), and with energies 0, oji and L02 respectively) 

H s =u 1 \l){l\s+u 2 \2){2\ s . (57) 

This system is also refered to as V-type three level system. This system is coupled to a bath of non-interacting harmonic 
oscillators at zero temperature: 

H B =J2^kb{b k . (58) 

k 

We choose the system-bath interaction such that transitions between the second and the third level are forbidden: 

fl>=|l)<0|Bi + |2)(0|B 2 +h.c, (59a) 
Bi=Y,9ikb k . (59b) 

k 

This is a generalization of the one-excitation example in Ref. fl]][p.461]. The "number operator" 

N=\l)(l\ + \2)(2\+Y,b{b k , (60) 

k 

commutes with the closed system Hamiltonian H = Hs + Hg + Hj and therefore is a conserved quantity. Let us define 
Hq = Hs + Hb and transform to the interaction picture: 

Hj(t) =e lHot H I e~ lHot (61a) 

=e Ml *|l)(0|B 1 (t) + e^*|2)(0|B 2 (t)+h.c. , (61b) 

where 

B, l (t)=Y,9ikb k e~ l ^ t . (62) 



As there are no excitations in the bath at zero temperature, we can restrict ourselves to the subspace with eigenvalues N = 0,1, 
which is spanned by \0) SB = |0) ® |0) B , |1) SB = |1) ® |0) B , \2) SB = |2) ® \0) B and \k) SB = |0) ® \k) B . Note that \k) B 
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corresponds to having a single excitation of the k-ih oscillator. Therefore, we can write the wavefunction at time t as 

^(i)) = co|0) + Cl (i)|l)+ C2 (t)|2) + ^; Cfe (t)|fc) , 



(63) 



where Co is time independent since H] (t)|0) = , Vi. For later comparison with the master equation results we need the reduced 
density matrix 



As = tr B (|V)<V'l) = 



1 - |ci| 2 - |c 2 | 2 c cl c c* 2 



C Ci 

c* c 2 



\Cl\* CiC 2 

c\c 2 |c 2 | 2 



where we omitted the time arguments for notational simplicity. The interaction picture Schrodinger equation 



(64) 



yields a set of coupled differential equations for the coefficients (time arguments omitted): 



, * — i(uj-\— cjz. )t , * —iiujo —uju)t n 



We can integrate Eq. d66bb and plug it into Eq. d66al which yields for c\ 



ci 



E 



(65) 

(66a) 
(66b) 

(67) 



Lets assume for simplicity that gu- = g 2 k = 3fe- Note that if lui = lu 2 , the state |— ) = -4= (|1) — |2)) is a dark state, i.e. it is 
decoupled from the bath, as both transitions — > 1 and — » 2 are coupled to the same bath operator B now, and the model 
reduces to the two level model from 1 1 1 . 

For the bath of harmonic oscillators, we will assume a continuum limit with spectral function J(lo) and a coupling constant g: 



(68) 



such that we can replace the sum over bath oscillators by an integral over the spectral function 

ci = -fi * c x - e^-^h * c 2 , 
C2 = -f-2*c 2 -e^-^ t f 1 *c 1 , 

where * denotes the convolution 



(69a) 
(69b) 



(/*<?)(*)= / dTf(t-r)g(r), 



(70) 



and 



f j (t) = e %u i t / dujJ(uo)e- wt = e w ' l f{t) 



(71) 



is the frequency shifted Fourier transform of the bath spectral function. Eqs. ( |69a| i and ( |69b[ ) can be solved numerically to give 
the exact solution for the time evolution of the system density matrix. We use a fixed step size fourth order Runge-Kutta method 
and calculate the convolution integral using the Riemann sum approximation with the stepsize of the Runge-Kutta method as 
interval width. 
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— L CO, t 
8000 



cj c t 



2000 4000 6000 8000 2000 4000 6000 8000 

(c) Coherence, \u)\ — u>2 \ = 10g (d) Coherence, \uj\ — u)2 \ = § 

FIG. 1: Exact dynamics of the populations and coherences for g = 0.001cj c and ps(0) = |1) (1|. In (a) and (c) we chose 
u>i = 0.095u; c and u-2 = 0.105w c , in (b) and (d) U\ = 0.09975w c and lj-2 = 0.10025cj c . In (a) and (b) black and dotted is poo, 
blue and solid is pu and red and dashed is p22- In (c) and (d) we display the coherence between levels 1 and 2, i.e. pi2- 



V. PARAMETERS OF THE MASTER EQUATIONS 



A. RWA-SME 



We now wish to derive the SME parameters associated with the model described in section|lV] First, we note that the interac- 
tion Hamiltonian in Eq. ( |59a[ ) is of the form: 

Hi = A ® B + A f ® B f , (72) 

where 

A=|l)(0| + |2)(0|, B = Y,9kb k - (73) 

k 

In the derivation of our SMEs it was assumed that A and B were Hermitian, but here they clearly are not. This does not change 
the overall form of Eqs. ( |48| l and ( |50"] l, but we must be careful how we define the bath correlation function. Since we assume 
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1000 2000 3000 4000 5000 6000 1000 2000 3000 4000 5000 6000 

(a) CG-SME (b) RWA-SME 



Pii 
1.0-tt 




1000 2000 3000 4000 5000 6000 
(c) Exact 

FIG. 2: Dynamics of the populations pa calculated with the different approximate equations compared to the exact dynamics at 
g = 0.001w c , LJi = 0.095u c , L02 = 0.105cj c . The dashed red line is p-nii), solid blue line is p\\(i) and the black dotted line is 
Poa{t). The optimal parameter for the averaged evolution generator is u c At = 63.7 for this set of system parameters, in 
agreement with Eq. ( fTo*| ). As initial condition we set ps{0) = |1) (1|. 



T = 0, the bath state is ps = |0) (0| and therefore we have the correlation functions (dropping the B subscript for simplicity): 

(Bt(t)B(0)) = (B\t)B\Q)) = (B(t)B(Q)) = (B\0)B(t)) = (B^0)B^(t)) = (B(0)B(t)) = , (74a) 

OO 

(B(t)£t(0)) =J29 k gl^{b k bl\0)(0\}e-^ t =^|p fc | 2 e-^* = / duJ(u)e** =: B(t,0) , (74b) 

k .k k q 

oo 

(B(0)B\t)) =Y / 9kgl,tr{b k b{,\0)(0\}e^ t = ^| 5fc |V<^ = J dw J{w)e iut = B(t, 0)* , (74c) 

k .k 1 k q 

where we assumed the spectrum of the bath to be lower bounded by 0. The one-sided Fourier transform of one of the nonzero 
correlation function is: 



oo oo oo oo 

T+H = J dtB(t,0)e tut = J dw'J(u>') J dte- u{uJ '- u) = J du'J(u') \ nS(u>' - u) + iV (jjzr^jj 



(75a) 



o n 

oo 

' 11. J\ 



=7rJ(o;)e(w) +i I dJ J{J)V ( — — - I , (75b) 
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where we used the fact that 



dte lu,t = h5{uj) + VP ( - ) . (76) 







where we have used the exponential integral function: 

oo 



(77) 



V denotes the Cauchy principal value and and denotes the Heaviside step function. Similarly, 

oo oo 

r_(w) = J dtB(t,0)*e tuJt =nJ{-uj)<d{-uj)+i J du'J(w')V ( _J l — J = T + (-u;)* . 
o o 

We assume our bath is characterized by an ohmic spectral density 

4 c (w) = r/uie~^ , (78) 
where r\ has the dimension of time squared. We can calculate the integral in Eq. |75b) to give: 

oo oo 

r) J dui'oj'e~^V ( - 1 ^ = - r]UJ c e~^ J du e~ u ^1 + ^-V QXj , (79a) 
o -g- 

= — r)Lu c + -quje u <= Ei I — , (79b) 



(80) 
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The function has branch points at z — and z = oo and the branch cut is defined to go from z = — oo to z = . Hence we get 
for the one-sided Fourier transform of the bath correlation function 



r +M = + iS(u) = 7rJ Wc (w) + i 



rju c + J UJc (u)Ei 



r_(-w) = -7(w) -iS(u) . 



Similarly, 



This allows us to write the generator from Eq. ( |48| ) as: 

= ~i [H L s,Ps(t)} + ^ 7(^j) ^ajp 5 (i)a] - ~ { at a j, 

where 



J'=l,2 



i7 L s = 51 S ( U3 j) a j a i ' 
i=i,2 



(81) 



(82) 



(83) 



(84) 



is the Lamb shift Hamiltonian and a, = |0)(i|. Transforming back to the Schrodinger picture just adds the system Hamiltonian 
to the unitary part. 



B. CG-SME 



Here we construct the parameters of the CG-SME, Eq. < |50] ). For the zero temperature case with the correlation function ( |74b| >, 
Eq. (Hob simplifies to (for detailed calculations see Appendix [P| 



J[y) 



*w(f) = I dv- "-rr-. r f e" v 

' \ui — v)[u' — v) \ 



i((jj' — oj)t —u)t — i(uj — f)t 



For lu 7^ <J this yields 

*w(t) = \(l + e *(«'-«-)*) (h(u) - h(u')) + e*"'* - J 2 ») + e"^ (J a (w') - J 2 (w)) 

— CJ L \ / 

where for the ohmic spectral density in Eq. ( |78] i: 

I\(l>j) = —goj c + guie^^Ei 



Uu) = ^— + gue-^-^Ei ( u (—-it 

1 — iuj c t 



For u = u', Eq. ((85|) reduces to 



/,_._.(/) = 2 / d 3^±^ (l - COs(l/*)) 



For the ohmic spectral density this yields 

+277(1 - cos(wt)). 
Knowing the form of & we can now compute the rates 7 in Eq. (|51a|i. 



(85) 



(86) 

(87a) 
(87b) 



(88) 



1 - — - iut ) Ei ( — + iut I + ( 1 - — + iut ) Ei ( — - iut ) + 2(— -l)Ei( — 

UJ r: \Ur \ Ur \U C . / \U n I \U C , 



(89) 
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0.02 





20 50 100 200 300 

(a) uj = 0.1w c , Suj = O.OluJc, g = O.OOkJc 
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(b) uj = O.lwc, Suj = O.OlaJc, 9 = 0.002cj c 
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(c) uj = O.lcJc, 5lj = 0.01u> c , g = 0.003w c 

D[ 0,2667] ( PAE,Pcx»ct) 

0.10 r 
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(d) w = OAuJc, Suj = O.OltJc, g = 0.001o; c 



Wr At 
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0.114 



0.02 



20 50 100 200 300 

(e) uj = 0.3uj c , Suj = 0.03o; c , g = 0.003w c 



6J. At 



FIG. 4: Integrated trace-norm distance of the solution of the CG-SME and the exact dynamics (green) for different averaging 
times At compared to the RWA-SME (black, dashed). The integration intervals are chosen according to the relaxation times 

that result from the different system parameters. 



The Lamb shift term for the zero temperature bath correlation function in Eq. (|74b[) is 



S U u'(t) = / dvJ{v) I ds 



s(cj'— uj) p is(oj — v) p is{u' — cj) p is(y—uj) 



iiv — uj) 



i(y — uj') 



(90) 



For uj 7^ uj' we get 

Suju)' if) = — ; 



(l - e «(«'-«)) (I.iuj) + h {fJ)) + e -** (J a (a/) - I 2 (uj)) + e-'* (7 2 » P 2 {uj')) 



2{uj-uj') 

where we have again used the integrals defined in Eq. ( [87) . For uj = uj' , Eq. (|90) reads 

Suj U (t) = -tli(u)+h(u) , 



(91) 



(92) 
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FIG. 5: Trace-norm distance of the different approximate solutions to the exact solution, for the same paramters as in Fig. [2] 
The dotted red graph corresponds to the RWA-SME solution, the solid blue one to the optimal Ai opt CG-SME solution, and the 
dashed green and dashed-dotted black line are for At = Ai opt /2 and At = 2Ai opt , respectively. 



where 1\ is defined in Eq. (|87a[) and 1% is 



J3M = v- 







) El \U C 


- iujtj 







2 1 



We now have everything we need to construct the CG-SME, Eq. d50 



2iojt I Ei ( iujt 



+ r]sin(ujt) . (93) 



VI. COMPARISON BETWEEN THE MASTER EQUATIONS AND THE EXACT DYNAMICS 



We first discuss the exact dynamics of the model presented in section IV For the calculations we set 77 = oj^ 2 . Fig.[T|shows 
two examples of the dynamics for different system parameters. Since [Hs, Hj] ^ for \uji — oj2 1 7^ oQthe dynamics are 
dissipative, so the system will eventually relax into the T = Gibbs state peabbs = |0) (0|- F° r l w i ~ ^2! ^ <?> the system relaxes 
into the T = Gibbs state very rapidly [see Fig.[TJa)]. For \u>i — u>2 \ ^9, two timescales can be observed: the initial population 
of |+) = -h= (|1) + 1 2)) decays rapidly, after that a strong oscillatory (with slow decay) behavior is observed in the |1) and |2) 
populations [see Fig. nib)]. 

We are interested in the regime where the RWA is no longer valid. Recall that the RWA is argued to hold when the time 
scale associated with |wj — w 2 | _1 is much less than the relaxation time scale ~ 1/g. Therefore, the RWA should start to 
break down when \ui — CJ2I > g- We investigate the behavior of our two SMEs in this regime. We choose the coarse-graining 
time scale A< such that it minimizes the integrated trace-norm distance between the CG-SME density matrix pcg and the exact 
dynamics pexact: 



£ , [0,t ma x](PCG, Pexact) 



1 



2/,« 



dt\\p C c(t) - Pexact(t)||l , 



(94) 



where || • ||i is the Schatten 1-norm (PHi := TiVAtA). In Fig. H we show an example of the dynamics of the different master 



1 Recall that as already pointed out in Section IV for u\ = ui2 we have a dark state that is decoupled from the bath. 



16 



3 r 
2 - 



0.2 - 
0.1 - 



1 



& LOCO ' i S UKi) ' 

3 r 
2 - 




i: ilA, 



0.5 

0.2 y 
0.1 - 



10 100 1000 

(a) ui = 0.1cj c , &ui = O.Oli^c 
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FIG. 6: The rate functions of the CG-SME. The dotted red line is b uiu>2 , the solid blue line b UlUl , the dashed-dotted red line is 
S UlUa and the dashed blue line corresponds to S U1U11 . The horizontal lines are the asymptotic value of the rate functions, i.e., 
the value they converge to for high oj c At. The vertical lines are at w c Ai opt . Note the qualitatively similar location of the 
w c Ai opt values relative to the features of the rate functions. This may be a clue for future studies in terms of independent 



determination of At 



opt- 



equations compared to the exact dynamics. In this example, we consider an initial state given by ps(0) = |1)(1|. The RWA- 
SME [Eq. |48|] shows no population in the state \2) since the Lindblad operators only allow for transitions between |1) — > |0) 
and 1 2) — > |0). This is in stark contrast to the CG-SME [Eq. (j50j] which properly shows the oscillatory population in the state 
1 2). We can compare the results quantitatively by computing the trace distances to the exact result, shown in Fig. [3] The two 
SMEs show similar qualitative time behavior in the trace distance, but the CG-SME outperforms the RWA-SME. We calculate 
the optimal At for the three level system for different sets of parameters (shown in Table|l|: 



VII. HOW CAN WE FIND A GOOD At FOR REAL PROBLEMS? 



If we want to make the CG-SME useful more generally, we need to find a way of determining At without comparing it to an 
exactly calculated result. We calculate the integrated trace-norm distance of the solutions of the CG-SME and the exact solution 
for suboptimal values adjacent to the optimum, shown in Fig.|4] 

The numerical result suggests that for this particular three-state model, the optimal averaging time At strongly depends on the 
timescale lo^ 1 of the transitions driven by the bath, and does not depend strongly on the other parameters g and 5lu. The latter 
are the ones that determine the validity of the RWA. On the other hand we can see, that if the RWA fails, i.e. if the RWA-SME 
solution has significant trace-norm distance from the exact solution, there is a wide range of values for the averaging time where 
the CG-SME does not yield optimum results but still outperforms the RWA-SME. This fact is further clarified in Fig. [5] where 
using half the optimal value as well as its double still outperform the RWA-SME on average. 

It is instructive to look at the rates of the CG-SME for different At (Fig. [6}. The optimal value of At is located at a 
similar point on (the similar looking) graphs, independent of the parameters of the model. This fact might make an analytical 
determination of the optimal averaging time possible, but we leave the investigation of such a possibility to future research. 
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TABLE I: Optimal At for different parameters. 
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VIII. CONCLUSIONS 



First-principles derivations of quantum Markovian master equations have so far relied on the rotating wave approximation. 
Here we have shown how the RWA can be avoided in a first-principle derivation, and replaced by a coarse-graining procedure. 
The coarse graining timescale is a free, phenomenological parameter which is a priori bounded between the bath correlation 
time (inverse of the high-frequency cut-off in the bath spectral density) and the characteristic timescale for significant changes 
in the system density matrix in the interaction picture, and can be optimized. To test the new, coarse-graining based SME, we 
applied it to a simple three-state model coupled to an infinite harmonic oscillator bath, whose exact dynamics we can solve for. 
Comparing, we found that the standard RWA-SME misses important qualitative features of the solution, while the CG-SME 
captures these features relatively well. By optimizing the coarse-graining timescale we achieved significantly better agreement 
with the exact solution by using the CG-SME rather than the RWA-SME. To perform the optimization we required the exact 
solution; however, we showed that there is a substantial range of the coarse-graining timescale parameter that still allows the 
CG-SME to outperform the RWA-SME. Future work should address how to extract reasonable and advantageous estimates of 
the coarse-graining timescale without relying on exact solutions. 
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Appendix A: From Eq. {37} to Eq. ((39) 



The unequal upper integration limits in Eq. ( |38a| > are inconvenient. To remove them we note the following relations for the 
integral, where for notational simplicity we suppress the i-dependence for now: 



pt pS |~ pt pt pt pt 

B^ = ds ds , e l ^' s ' uls ">B(s, s') , = / ds ds' - ds ds' 
Jo Jo Uo Jo Jo Js 



i{uj r s~ujs') 



ds I ds' — I ds' I ds 



B(s,s') 



e i( w ' a -«.') J B( ai s >) , 



(Ala) 
(Alb) 



= [ ds [ ds'e l{ul ' s - UJS ' ) B(s,s')- [ ds f ds'e l ^' s '-^B{s' » , (Ale) 
Jo Jo Jo Jo 



where & WL y [Eq. (|40j>] has the desired equal upper integration limits. It follows immediately that 



(Aid) 



(A2) 
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and from Eq. d40| we also directly get 



(A3) 



Therefore, Eq. p7] > yields terms such as: 



(A4a) 



(A4b) 



where in the second term on the RHS we have switched w a/. Furthermore, this second term is now exactly in the form of 
the original term, so we have the result: 



53 [A>pAl,B uu > + AspAl (B^)*] = o [ A "P A l> b ""' + A ^P A l h ^J\ = 53 A "P A l> b ^ 
For the remaining terms, we write: 

+ \Y [ A l' A "P ( B ^' - (B U ' W )*) + MU,' (ft/a, - (Bo*;')* 

= \ 53 + ~ o^r) + Ml- a 

Putting all this together we can now go from Eq. ([37} to Eq. < (39] >. 



(A5) 



(A6a) 
(A6b) 



Appendix B: Relation between the Coarse-Grained Master Equation and the Standard Lindblad Equation 

In order to study the relationship between the CG-SME and the RWA-SME, we begin with the dissipative coefficients of the 
CP map: 



«w(i) = da ds'e 1 ^' B(s,s') . 



We switch to new variables u = s — s' and v = s + s' using the homogeneity of the bath correlation function: 

t v 2t 2t-v 

bu U >(t) = ^J dve'^ v J due i ^ u B(u,0) + dve l ^ v J due l ^ u B(u,0) . 

-ii t 



-(2t-v) 



We replace v by 2t — v in the latter integral: 



W(*) = \ f dve l ^ v ( due l ^ u B{u,0) + ^ 



dye' 1 — v / due 1 — u B{u,0) 



(Bl) 



(B2) 



(B3) 



(I reus I - - U (v-t)) I due l ^ u B(u,0) 



(B4) 
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For lu = <J this expression reads 



h^(t) = J dv J due wu B(u,0) , 

-v 

and partial integration yields 

t 

b uu (t)=- [ dvv (e luJV B(v,0) + e- %UJV B{-v,0)) +t I due iuu B(u,0) 



= t dull- 
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e^ u B(u,0) . 



(B5) 



(B6a) 



(B6b) 



For our SME, we are interested in b uu r(At)/At, Setting At — > oo, the part proportional to £fB(u, 0) is negligible, and we 
recover the RWA result of Eq. d49l: 



lim ^-b uu (At) = I due l " u B(u,0)=~f(u) 

At— S-OO LAI 



(B7) 



For ui 7^ <J we also perform a partial integration, but this time the boundary terms vanish: 



' (a/ - uj) 



, (uj 1 — uj) . . 
dvsin I -(v - t) 



e l ^ v B(v, 0) + e~ i ^ v B(-v, 0) 



(B8) 



Changing from v to — v in the second term we get 



, , . 2e*T* 

*W it) = - — r 

(uj' — Uj) 



dv sin 



(w' - u) 



(v-t)\ e z — v B{v,0) 



o 



dv sin 



e 2 



(a/ - u) 



{-v-t) e z — v B(v,0) 



dv 



sm 



t I ( e 1 ^ + e lu v 



sgn(v) / u)' — ui 



cos 



(e iuv 



(B9a) 



B(v,0), (B9b) 



where we used the angle sum identity for the sine in the last equality. If we again look at the term b^i (At) j At and take the 
limit of At — > oo, we get 



lim : = lim 



At->oo At Ai^oo (ul 1 — U})At 



sm [ -At ) (7(0;) +7(0;')) + 2 cos 



At (S(uj) - S(oj')) 



(BIO) 



where we have used that for T(u) = J °° dse luls B(s, 0), we have 7 (w) = r(w) + T(uj)* and 2iS(w) = - T(ui)*. Since 
nothing cancels with the overall At -1 , we find that the uj ^ uj' term vanishes, which is what we would expect for the RWA 
result. A similar calculations could be done for the Lamb shift term ( |42| >. Therefore, the RWA results can be understood as the 
At — > 00 limit of the coarse-graining parameter. 
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Appendix C: Calculation for the Two Level Example 



by: 



We consider the model in section III Let us denote the eigenvalues of H$ by e± = ±wo/2 and their respective eigenvectors 

(CI) 



|e+)^(l,0) T , |e_)^(0,l) T , 



and T denotes the transpose. Note that the Lindblad operators in both the CG-SME and RWA-SME are given by: 

A- Uo = |e+)(e+|CT 2 |e_)(e_| =0 

A Q = \e + )(e + \a z \e + )(e+\ + \e^)(e^\a z \e^)(e^\ = a z 



| e _)<e> x | e+ )< e+ |=0 



(C2) 
(C3) 
(C4) 



so in this two state model, only the uj = term contributes to the sums over u. This in turn gives that H^s oc 1 such that 
[-HlSj Ps(t)] = 0. In order to solve the SMEs, it is convenient to vectorize the master equation. Let us define: 



p = vec(p) 



Pi2 

V i J 



where vec(p) corresponds to stacking the columns of p (in some basis). We now use the identity 1 16 1: 

vec (ABC) = (C T ® A) vec (B) 
where (A, B, C) are arbitrary matrices, such that we can write both SME as 



vec ( —ps ) = 7 ( A* ® A 



t T ® A\A 



o + (4^o) T ® tYj Ps = Cp s {t) 



where it is to be understood that 7 = 7cg for the CG-SME and 7 = 7rwa for the RWA-SME. These are defined as: 

7CG = At J d ' S J ds ' i B ( s > s ') + B ( s 'i s )) ' 7rwa = y dsB(s,0) 
The solution written in terms of the superoperator C is simply 

p s (t) = exp(£t)ps(0) 

Using that A a = a z , C reduces to : 

C = 7(0-3 ® a z -1 ®1) 
Conveniently, C is diagonal with entries (0, —27, —27, 0), so we can immediately write: 

(l \ 

exp(£) = 



exp(— 2ft) 



exp(— 2ft) 



V 



1/ 



Therefore, we find as a final result that: 

Ps(t) 



pn(0) exp(-2 7 t) / 5 12 (0) 
exp(-2 7 t)p 21 (0) p 22 (0) 



(C5) 



(C6) 



(C7) 



(C8) 



(C9) 



(C10) 



(Cll) 



(C12) 
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Moving to the Schrodinger picture, the result is adjusted to 

p„( t )=( P 11 ^ exp(-2jt - iu) t)pi 2 (0) \ 

\ v exp(-2 7 t + ^)P2i(0) p 22 (0) )' 

Finally, for the infinite harmonic oscillator bath to bath operator B = X ^ k {bk + bl), we can explicitly calculate the form of the 
7's by using the result for the two point correlation: 



(B(s)B(s')) = (B(s ~ s')-B(O)) = A 2 J du - J _^l fju) (e^ 8 "^ + eM'-^-P^ 

This then gives: 

2 1 f°° T/ , , fl -cos(cjAt)\ - f 00 . , , //3uA Aisinc 2 (^) 



(C14) 



^ G=2X AtJ ^)= X ^J J (-)oot^) ^ (C15a) 

7 rwa =2ttX 2 lim ( ^ (C 15b) 



w->o \1 - 

When 7cg is written in the form of the second equality, it becomes clear that taking the At — >• oo limit corresponds to the RWA 
result since we can use the identity: 

sine 2 (^) 

lim At- \^-L = 5(u) (C16) 

At— >oo 7T 

and 7cg becomes 7rwa- 



Appendix D: Calculations for the Three Level Example for the Averaged Evolution Generator 

Here we provide the detailed derivations of the expressions for the averaged evolution generator of the three level model in 
section [IV] Let us calculate Eq. ( |40] > for the zero temperature case with the correlation function Eq. ( |74b| >: 



W(*) = ds ds'e 1 ^' 3 -^^ / dve- iy{s ~ s,) J(v) 



dvJ{v) \ dse ls ^'-^ / ds'e- is '^- v ^ 



dv 



J(y) 



(u — v)[bj' — v) 



(j(u'-v)t _ f e -i(w-v)t _ 



dv 



(u — v)(uj' — v) 



(Dl) 



For the second equality sign we used the fact that the Fourier transform of the spectral function converges absolutely for a well 
behaved bath: 



dv 



oo 

e -iu( a -,')j^ = J dv \J(y)\ <o0 




(D2) 



Let us look at the case u> ^ io' . With the help of 

1 



(u> — v){lu' — v) id' — LJ \uj — v oj' — V ) ' 



(D3) 
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we can split up the integral: 



?w(f) 



Of — UJ 



dvj(y) 



w — f or — f w — v uj' — v 



OJ' — V UJ — V UJ — V 



(l + e^-**) (h(u>) - h(uj')) + e-'* (^(o/) - 7|M) + e"*-' (J 2 (o/) - J 3 (w))' 



where we have defined the integrals 



(D4) 



(D5) 



For the ohmic spectral density d78) the integrals 1\ and I2 are analytically solvable: 



Ii{uj) = rj / dv 



rj j dv 



(v + uj)e 



-rjuj c + rjuje <*c Ei 



I<z (oj) = rj I dv = —rj I dv 



(D6) 



■ rjuje v ) 



1 — 



Ei uj 



it 



(D7) 



where we substituted v by + uj in both calculations and used the definition ( |80) . For the case uj = uj' we get from ([85}, again 
substituting vby v + uj 



b uu {t) = 2 / dv . Jiy \j l ~ cos((^ - v)t)) = 2 / dv Jiy + "\ l - cos(i/t)) 





For the ohmic spectral function this yields 



(D8) 



b uu (t) = 2 V I dv ^ V + U) \ "° (1 - cos(z^)) 



rye 



Ei 



icji + Ei iuot — 2Ei 



2rjuje~~^ / di/ — ,5- (1 - cosfz/i)) . (D9) 



Partial integration simplifies the last integral to 



dv " (1 — cos(vt)) = I dv ( tsin(vt) (1 — cos(^t)) I + - — (1 — cos(cji)) 



— Ei 

UJr \ UJ C 



1 / 1 



2\uj, 



1/1 



it ) Ei ( h iut ) - - I it ) Ei ( iujt ) H (1 - cos(wi)) , (D10) 



so that we can write altogether 



b uu (t) = rje »■ 

+277(1 - cos(wt)) 







) E1 (w c 


- iujtj 







1 h iujt Ei iujt 



1 Ei 



(Dll) 



Let us now calculate the Lamb shift term for the zero temperature bath correlation function Eq. ( |74b| i. Similar to the case of 
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S UU '(t) = -\\ dvJ{v) I ds ( ds' ^ e *(«'»-««') e -<«'(«-*') _ e -i(u S -w> s ') e -iv(s'- s )^ 





t 



dvJ(v) / ds 



is{uj'-u) I £ g i gis' '(v-w) _ e is{v-u) / e is' {u' -u) 



dvJ(v) / ds 



dvJ{y) \ ds 



i I ns i - ( •• I as e 

o o 

,ls(u/-I>)~ * _|_ gls(!/-Oj)^ ^ 

— w) — a/) 



(D12) 



o o 

Now we have to distinguish the two cases, u = lj' and cj 7^ a/. For cj 7^ u/ we get 

O it(o/— u) _ ^ g— it(r/— ui') _ ^ git(u'— w) _ ^ ui) _ J 



Su,u>(t) = -- / dvJ(v) 



dvj(v) 



(lj — uj')(v — uj) (y — w){y — u') (u — u'){y — u)') (v — oS){y — w') 



(lj — uj')(v — uj) (lj — lj')(v — lj') (y — (jj)(v — lj') 



With the help of ( |D3[ ) and ( |D5| l this formula simplifies to 



2(w-u') 
For the case w = lj' Eq. d90| reads 



1 - e u ^'-"A (h(u>) + h{Lj')) + e-^ (I 2 (u/) - J a (a;)) + e^' 4 (I 2 » - 7|(a/))' 



(D13) 



(D14) 



^(t) =-7, dvJ(y) \ ds 



— uj) i(y — lj) 



dvJ(v) 



2t e i(v-u)t _ e -i(v-iS)t 

i(uj — v) (uj — v) 2 



dvJ(v) 



t sin((^ — uj)t) 



v — uj (y — lj) 2 



= -th(Lj) + I 3 (UJ) , 

where Ii is defined in ( |D5| > and I 3 is defined as 

00 

I 3 (uj) = [ dv,](v) 



sin((i/ — uj)t) 
(v-uj) 2 



(D15) 



(D16) 



For the Ohmic spectral density this yields 



. sm((v — oj)t) f , , . _j^sin(ut) 

7 3 (w) = 77 / awe »« — ? — — = rje w = / d-u(u + a;)e «c — 

(v — lj) z J u z 



. Ei I h iujt ) — Ei [ icjt ] ] +uie "c IUoj) 

2l V \LJr J \UJr. ' 



(D17) 
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where we defined 



Ii(ui) = I du 



e <*c 



s'm(ut) 



(D18) 



Partial integration yields 



= — / du 



du- 



e »< sin(wt) + t cos(ut)^J 



e^c sin(wi) 



<jj 



e "° 


e itu , 


'jL . 




















2u c 





sin(wi) 



2 //. 



1 \ , f uj \ ( 1 \ . / a; \ sin(W) 
H Ei iwi + I i Ei \-iojt) -\ — - 

lW P J\U)r. I \ 10J C ■ \Wc J iO 



Using this result, we get the final result for /3 



J 3 (w) = ?y- 



2/ 



I— 

UJ r 







MS- 






- iuiij 



I 



0J r 



(D19) 



?7sin(u;i) . (D20) 



[1] H. Breuer and E Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2002). 

[2] C. A. Rodriguez-Rosario, K. Modi, A.-M. Kuah, A. Shaji, and E. C. G. Sudarshan, Journal of Physics A: Mathematical and Theoretical, 
|41, 205301 (2008)1 

[3] A. Shabani and D. A. Lidar, |Phys. Rev. Lett., 102, 100402 (2009)] 

[4] W. T. Pollard, A. K. Felts, and R. A. Friesner, "The Redfield Equation in Condensed-Phase Quantum Dynamics," in Advances in 
Chemical Physics (John Wiley & Sons, Inc., 2007) pp. 77-134. 

[5] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom-photon interactions: basic processes and applications Wiley-Interscience 
publication (J. Wiley, 1992). 

[6] C. Gardiner and P. Zoller, Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applica- 
tions to Quantum Optics, Springer Series in Synergetics (Springer, 2004). 

[7] D. Walls and G. Milburn, Quantum Optics SpringerLink: Springer e-Books (Springer, 2008). 

[8] P. W. Milonni, J. R. Ackerhal t, and H. W. Galbraith,|Phys. Rev. Lett., 51, 1108 (1983)| 

[9] K. Zaheer a nd M. S. Zu bairy, Phy s. Rev. A, 37, 1628 (1 988) 
[10] M. D. Crisp, |Phys. Rev. A, 43, 2430 (1991)| 

[11] G. Ford and R. O'Connell, Physica A, 243, 377 (1997) 

[12] C. Fleming, N. I. Cummings, C. Anastopoulos, and B. L. Hu,|J. Phys. A, 43, 405304 (2010)1 
[13] D. Bacon, D. A. Lidar, and K. B. Whal ey,|Phys. Rev. A, 60, 194 4 (1999)| 
[14] D. A. Lidar, Z. Bihary, and K. Whaley, Ch em. Phys., 268, 35 (2001)T 
[15] R. Alicki, [PhysTRev. A, 40, 4077 (1989)| 

[16] R. Horn and C. Johnson, Topics in Matrix Analysis, Topics in Matrix Analysis (Cambridge University Press, 1994). 



